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We analyse the dynamics of trapped matter shells in spherically symmetric inhomogeneous A- 
CDM models. The investigation uses a Generalised Lemaitre-Tolman-Bondi description with initial 
conditions subject to the constraints of having spatially asymptotic cosmological expansion, initial 
Hubble-type flow and a regular initial density distribution. We discuss the effects of shell crossing 
and use a qualitative description of the local trapped matter shells to explore global properties of 
the models. Once shell crossing occurs, we find a splitting of the global shells separating expansion 
from collapse into, at most, two global shells: an inner and an outer limit trapped matter shell. 
In the case of expanding models, the outer limit trapped matter shell necessarily exists. We also 
study the role of shear in this process, compare our analysis with the Newtonian framework and 
give concrete examples using density profile models of structure formation in cosmology. 
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I. INTRODUCTION 

Studies of non-linear structure formation in cosmology, 
namely spherical top hat collapse models, often assume 
that there is no influence of the cosmological background 
on a finite domain which has disconnected from the back- 
ground dynamics (see e.g. 

We have looked at this problem in more detail in 
Ref. Q and found local conditions under which such sep- 
aration could be justified for inhomogeneous cosmologi- 
cal models. In particular, we have studied the possibil- 
ity for perfect fluid solutions to exhibit locally defined 
separating shells between collapsing and the expanding 
(cosmological) regions. 

The simplest examples given in [f| were set in an inho- 
mogeneous universe of dust with a positive cosmological 
constant and the nature of the dust spherical shells al- 
lowed the system to be entirely determined from its initial 
conditions, at least, until the eventual occurrence of shell 
crossing. 
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However, shell crossings or caustics are expected to 
happen in these settings with more general initial condi- 
tions than in [|| and an interesting question is whether 
our previous results are robust with respect to the occur- 
rence of shell crossing. This is the main concern in this 
paper, which can be regarded a natural follow up of our 
previous work Q . 

Shell crossing in spherical symmetry has already been 
studied in several past works, although in contexts differ- 
ent from the one of present paper. For Lemaitre-Tolman- 
Bondi (hereafter LTB) spacetimes shell crossing condi- 
tions were established by Hellaby and Lake Q in terms 
of the metric data and more recently re- written by Suss- 
man in terms of quasi-local scalars H Q . Goncalves , 
has shown that shell crossing exists for A-LTB spacetimes 
with charge. In |Tlj], it has been shown that shell-crossing 
occurs for a large class of initial conditions in models of 
formation of voids and some cases of fluids with pressure 
gradients. 

There were also several works about the strength of 
shell crossing singularities, with the general conclusion 
that it is a weak singularity in the sense of Tipler [l^ . 
This then raised the question of the continuity of the 
metric across these singularities and, very interestingly, 
solutions of dynamical extension through shell crossing 
sing ularities of LTB have been proved to exist, by Nolan 
|13| , while the case including a cosmological constant and 
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electric charges has been discussed by Gongalves [To| . 

A complementary treatment was given by Nunez et al. 
fl4j for metric extensions through shell crossing based 
on the interactions between shells, which translate in a 
conservation relation between mass and momenta, for 
timelike massive shells. Physically, this conservation re- 
lation summarises the microphysics of the fluid, however 
for dust, only purely gravitational interaction occurs be- 
tween cros sing shells, hence the rest mass of each shell is 
conserved |15| . 

Here, we shall not deal with the problem of metric ex- 
tensions after shell crossing and, motivated by the above 
results, we shall assume the validity of the field equations 
in between shell crossing events and the continuity of the 
radial coordinates. Our main concern here will be to 
study the effects of shell crossing on the existence and sta- 
bility of separating shells in spherical symmetry. In this 
paper, we shall also discuss the role of shear in the forma- 
tion of shells which separate expanding from collapsing 
regions, we shall compare our results with Newtonian 
cases and give a concrete example of initial data which 
develops shell crossing and exhibits separating shells in a 
A-dust model. 

The models considered in this paper obey the following 
properties: (a) spherically symmetric dust (the rest mass 
of infinitesimal pressureless shells is conserved under shell 
crossing) with a cosmological constant in Generalised 
LTB (GLTB) system ; (b) Lagrangian treatment of the 
radial coordinates (assume there are metric extensions 
through shell crossings); (c) asymptotic spatial cosmolog- 
ical behaviour (Friedmann-Lemaitre-Robertson- Walker, 
hereafter FLRW, at spatial infinity); (d) initial Hubble- 
type flow (outgoing initial velocities); (e) regular initial 
density distribution (no finite mass for infinitely thin 
shell, and no singularity or zero density at the centre). 

The paper is organized as follows: in a first part ((III) we 
recall the conditions for the existence of matter trapped 
shells and study the role of shear on the existence of those 
shells in A-LTB models. Section (jllip is devoted to the 
study of the effect of shell crossing in A-LTB models. In 
particular, we perform a dynamical analysis and separate 
this study into a local and global effects. We give con- 
crete examples in section (|IV[) before presenting the final 
conclusions. 



II. TRAPPED MATTER SHELLS IN A-CDM 



A. Conditions for the existence of trapped matter 
shells 



In this section, we briefly recall some results of our 
previous paper Q which did not consider shell-crossings. 

The GLTB system proposed in Refs. @,[H|, has the fol- 
lowing simple form for the case of a A-dust model where 
P' = and P = P dust = (here we set G = 1 = c, A > 0, 
a is the lapse function, r (T, R) the areal radius and E 



the energy of spatial hypersurfacetQ) 

2 - + (1) 



ds 2 = -a(t,Rydr 



The Bianchi identities projected along and orthogonal to 
the timelike flow n = dt yield (P is the pressure, p the 
density, the prime / denotes Or, a dot ' stands for dt and 
O is the expansion along the flow) 

P = -(p + p) 3 e, (2) 

P' a' 

=— = => udt = dt* => a = 1, (3) 

p + P a > W 

and the Einstein Field Equations (M is the Misner-Sharp 
mass defined as M = ^ 47rpr 2 r'dR) 



Er ' = 2ri±|p' = T 2— F 'm/2- + -Ar 2 + E 
p + P p + P V r 3 

(4) 

=> E =0, E = E(R), unless there is shell crossing, (5) 



M = — MnPr 2 = +4nPr 2 a\l 2— + \kr 2 +E = Q 



r 3 



M =M(R), unless there is shell crossing, 

r 2 =2 h -Ar 2 + E. 

r 3 



(6) 
(7) 



(8) 



Efr ime derivation of Eq. ([5]) gives a Raychaudhuri equa- 
tion related to the gTOV function of Ref. [f| : 



gTOV =% - = -f. 



(9) 



The dynamical analysis detailed in Ref. [5j] yields the 
motion of separated non-crossing shells in their respective 
effective potential 



2M A , 

E=V(r) = -r 2 , 

r 3 



(10) 



where the unstable saddle pointjfor which gTOV = 0, 
gives a local separating shell (see figures 1 and 2, and 
repeated in Fig. [1]), in the case when the shell's energy 
reaches its critical value. This separating (or "cracking", 
by analogy with Herrera et al. [18j|) shell is characterised 
by 



ru 



3AI 
~A~' 



Slim = -(3M)3 A3 



-Ar 2 



(11) 
(12) 



1 Actually, 3 R = —2 ^ E J^ so E is related to the 3-curvature. 

2 In case of shell crossing, E can be nonzero as r' = and M gets 
changed by the loss or gain of the mass from infinitesimal shell 
crossings, so E = E(t, R) and M = M(t, R), in that case. 
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Figure 1. Effective potential kinematic analysis (left) and phase space analysis (right) from pj. The kinematic analysis for 
a given shell of constant M and E depict the fate of the shell, depending on E relative to Eum. It either remains bound 
(E < < i?ii m ) or escapes and cosmologically expands (_B> > Eu m )- There exists a critical behavior where the shell will forever 
expand, but within a finite, bound radius {E = En m , r < rn m ). The maximum occurs at ru m = ^/3A//A. The corresponding 
phase space behaviour follows, the scales are set by the value of ru m = ^/3M/A while the actual kinematic of the shell is given 
by E. 



while the energy follows 

E =r 2 + V(r). 



(13) 



Definition 1. Local trapped matter shells in A-LTB are 
defined in GLTB coordinates as the locus i?* such that 



e 

3 +a 



^ = 
r 



0. (14) 



This definition follows from Eqs. (3.11) and (3.16) of [j| 
applied to dust with A. 

In A-LTB, conditions (|14p are reached by shells 
at time-infinity which are characterised by Eqs. ([8]) 
and (fT2"j) so that (see footnote jp E(t = oo,R+) = 
Eu m (t = oo, i?*) (defining R+) i.e0 



=2 



r 3 (T,JJ*) 



1 (3M(g^Aa 
3 A_ r^TjQ ■ (15) 



So, since here the Misner-Sharp mass M and energy E of 
each shell is conserved in time (without shell crossing), 
and E is thus set by initial M(R) and fi(R) profiles, one 
can characterise local trapped matter, or limit, shell by 
the intersections E = Eu m (see [f|, for details). Global 
shells emerge from the neighbourhood behaviour around 
those intersections which local study we give in Sees. MI XI 
and HITCl 



Erratum: Eq. (3.14) of has a sign typo. It should read 



gTOV 



(-) 



+ a 



Before studying the occurrence of shell-crossing we will 
now examine more carefully the role of shear in these 
settings. 



B. The role of shear in the existence of trapped 
matter shells 



In Ref. Q, we derived the relation between expansion 
and shear (see Eq. III. 10) and found that, in the presently 
studied model, the shear could be put in the form 



r 3 



2E 
r 



r' 



3M 
r 



(16) 



In the latter equation the terms within the brackets mea- 
sure the departures from the profiles E = E(t) r 2 and 
M = M(t) r 3 that one would expect from a homoge- 
neous, uniformly curved models . Indeed, in FLRW mod- 
els E = kr 2 and M <x p{t)r 3 . Moreover, we should stress 
that Eq. (fTB")) yields the shear in terms of non-local (in- 
tegral) quantities (E and M). We can now evaluate the 
expansion and shear at the limit shell defined by setting 
Eqs. §8§ and © to zero at time infinity in Ref. [5j|, which, 
with the conservation of E and M, is simply defined by 
Eqs. and (fT"2")) . Combining those equations yields 



E' 



2M' 



(17) 



First on the limit shell we can write, setting E = En m , 
\ 2 ML (l _ .Jj) + (i _ m^)\ 

i r' \r Tn m J r V r / J 



a = T- 



2 ^+r*-3r* m 



(18) 



With the definition of mass issued from Ref. [3'sEq. 11.27 
in GLTB coordinates so 

M' =A<Kpr 2 r', (19) 

we then express the shear of the limit shell as 



an 



=F 



\ 



A 



"1 - 



A I n 



(20) 



which in the limit of time infinity simplifies into 

A-47rp(r iim ) 



ft// 



3VA~ 



(21) 



This quantity does not vanish in general. Since at that 
locus we have = 3 — a), the expansion then reads 



e, 




which in the limit of time infinity simplifies into 

A - 4irp(r Hm ) 



±- 



Va 



(23) 



We shall now use a particular form of initial data in order 
to study in more detail the role of shear in the appearance 
of the diving shell. In the examples below we shall assume 
M > 0, p > 0, A > and E < around the origin. 

So, consider analytic initial data for A-LTB as in (l9l . 
I2f 



M(R) = R 3 miR\ m > 

i=0 
oo 

E(R) = R 2 J2E l R\ E o <0 



(24) 



then, from the expressions above, we derive 



ft// 



l (R) = ±A 1/ 



1/2 I 1 2 I I 2m l d , / 3m 



R+(- 



•3 Q2/3 I 2/3"" 1 v 2/3 



1 )i? 2 + 0(i? 3 



5/3 



E Hm (R) = -3§A3 [ m fiZ 2 + J^£ 3 + 0(.R 4 ) j (25) 

V 3?77, n / 



also, for the re-scaling r(to,R) = R, we get an expression 
for the initial shear distribution as (see also (2l|): 



£i + 2mi 



± - 



6A 

l^ 2 + 4m 2 (E 1+ 2 mi f\ R2 + om 



6 V -4 



2A 3 



with A(R) = ^£ , + A/3 + 2m . 

It is interesting to see that for a fixed shell R near 
the centre, bigger M (i.e. bigger 771,3) means smaller 
initial shear but bigger |i?/i m | and ru m for that shell. 
On the other hand, since bigger initial shear implies 
smaller |-E/i m | (i.e. smaller departures from Eu m = 0) 
and smaller rn m , one can argue that, at least around the 
origin (and for the above initial data), shear contributes 
to the appearance of "cracking" limit shells. This is in 
agreement with the results of Herrera et al. [l8j]. We 
summarize this result as: 



Result 1. Consider a neighbourhood U of the origin 
where the A-LTB initial data can be written as (|24p . 
Then, bigger values of the initial shear \a(to, R)\ in U, im- 
ply smaller |i?ji m | and favour the occurrence of trapped 
matter shells in U. 

For data which is asymptotically Friedmann at infinity 
we take functions which, at infinity, can be expanded in 
the forrrB 



+00 



M(R) = ^2rriiRT, E{R) =^EiR- 

i=l i=l 

with mi 7^ and E\ 7^ O.By taking asymptotic expan- 



i=0 



4 This data ensures that the solution approaches FLRW at the 
origin which is therefore regular. 



5 Note that we only assume this data form at infinity and not 
around the origin. Otherwise, we would have a non-regular ori- 
gin. 
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Figure 2. Overcoming local configuration of E intersecting 
Eiim. Phase space and effective potential trajectories from 
dynamical analysis of [j| give the local qualitative behaviour, 
emphasised on the radial axis. Inner shells on bound tra- 
jectories and outer shells on unbound paths forecast no shell 
crossing locally. Considering Eu m as corresponding to the 
Newtonian zero radial velocity axis in [24Ll25l|. this configura- 
tion is analogous to, e.g., figure 1 of [2q |. 



sions we find : 



if) K* + 5(sJ +0 <s> 



3ml 



t n2 , 2 m 2 r>i , 2 m 3 ~ < 1 



Ei im (R) = -3' m »R> + -^-R* + ^^ + 0{— ] 
\ 6 ml 6 ml R 4 

(26) 



while the initial shear is: 



a{t Q ,R) = =F 



E * 1 I Q( \ 

2x/3\/3£i + A + 6mi i? i?f 



So, again, bigger values of the initial shear \a(to, R) \ near 
infinity, imply smaller |-Eh to | and favour the occurrence 
of trapped matter shells. 

We shall return to this issue in section IIVI where we 
study other examples in more detail. 



and the density diverge^. In geometrical terms, at shell- 
crossing there is a discontinuity both in the extrinsic cur- 
vature Kij and in the spacetime metric. For the space- 
times considered here, those discontinuities are finite and 
the magnitude of the jump in Kij can be read from the 
expressions derived in 1^, [23[ . 

Hellaby and Lake [6| (see also 0|) have derived nec- 
essary and sufficient conditions for the occurrence of 
shell-crossing in LTB, in terms of the free initial data. 
Other works have used other type of conditions which 
are sufficient to avoid shell crossing an< i therefore imply 
9rt 7^ 0. For example, in the case of LTB, Landau and 
Lifshitz [llj simply assume ObT > and, in [28[, Hellaby 
and Lake impose the condition for a simultaneous big 
bang in their local analysis around the initial singularity. 

Here, we shall take a different point of view and write 
sufficient conditions for the occurrence of shell-crossing 
in terms of the local behaviour of M and E in the neigh- 
bourhood of some intersection, when it exists, of the en- 
ergy E with the critical energy Eu m . In order to do 
that we first observe that two local configurations are 
possible in the neighbourhood of the intersection: either 
E'>E' Hm orE'<E' Hm . 

In the case E' > E' Um , shells just inside the inter- 
section radius will have a lower E than their respective 
En m and will therefore be trapped in closed trajecto- 
ries, following the dynamical analysis presented in Fig. 
Q] In that case, shells just outside the intersection will 
display higher E than their respective Eu m and will ac- 
cordingly be free to escape to infinity on unbound trajec- 
tories. That shell distribution will lead to the separation 
of neighbouring shells, those inside the intersection being 
bound to a finite region while those outside will escape 
to infinity. This case doesn't entail neighbouring shell 
crossings and is presented on Fig. [5] 

On the contrary, in the case E 1 < EL m , shells just 
inside the intersection will have a higher E than their 
respective Eu m and will accordingly be free to escape 
to infinity on unbound trajectories whereas shells just 
outside the intersection will display a lower E than their 
respective Eu m and will therefore be trapped in closed 
trajectories. Because of the configuration of that shell 
distribution, shell crossings of neighbouring shells occur: 
those inside the intersection escaping to infinity will have 
to cross those outside which are bound to a finite region. 
This case is presented on Fig. [31 We summarize this 
result as: 



III. SHELL-CROSSING AND TRAPPED 
MATTER SHELLS 



Result 2. Let A 



Ei, 



and consider a A-LTB 



spacetime where there is i?* such that Aj/^ =0. Then, 



A. Sufficient conditions for shell-crossing 

In terms of the comoving coordinates of metric ([lj, 
shell-crossing is defined as a surface for which dnr = 



6 There can exist cases where djir = and the density does not 
diverge. At those regular extrema, the extrinsic curvature is 
discontinuous while the metric is continuous and finite (6, 22]. 

7 A comment on the occurrence of caustics when using a syn- 
chronous reference frame can be found in Ref. |27H (§97). We 
must notice though that the latter assumes that the strong en- 
ergy condition holds, whereas in our present case the cosmologi- 
cal constant evades that assumption. 
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Figure 3. Undercoming local configuration of E intersecting 
Eu m . Phase space and effective potential trajectories from 
dynamical analysis of [j| give the local qualitative behaviour, 
emphasised on the radial axis. Outer shells on bound tra- 
jectories and inner shells on unbound paths will lead to local 
shell crossing. Considering Eum as corresponding to the New- 
tonian zero radial velocity axis in [24l . [25fl , this configuration 
is similar to, e.g., figure 2b of [2^|. 



a sufficient condition for the existence of shell crossing is 



We point 
to E' < 
ered in If 



out that, for A = 0, our condition leads 
which is the condition implicitly consid- 
In that case, we simply obtain E[ im = 

= 0. In this sense, our sufficient con- 



-SrrprV 

dition generalises, for A / 0, the result of Ref. p, [Dp. 
We also note that their condition on bang times tb(R) is 
t'b(R) — 0j while we can also allow for t' b (R) > as long 
as tb(R) is less than the initial time to considered here. 

There is an interesting analogy between our shell- 
crossing condition and a similar condition in Newtonian 
theory. In fact, the Newtonian approach used in [U, [25j 
considers kinematic configurations in velocity-radius two 
dimensional phase space which lead to one (and three) 
dimensional Zel'Dovich pancakes (see Refs. |25l [29|, |30j 
for the classical cosmological spherical context). Their 
behaviour is similar to the local evolutions of the dynam- 
ical configurations in Figs. [2] and [3] While in 0, HE] the 



8 For LTB with A = 0, we recall that the necessary and sufficient 
conditions for no-shell crossing in [(| are: 

T' B < 0, E' > 0, M' > 

where Tb is the bang time, while the necessary and sufficient 
conditons for no-shell crossing in [?| are: 

T' B < 0, E' > 0, M' > 0. 

Therefore, if one of these conditions fails then there will be shell 
crossing. 

9 Although our interest is in the neighbourhood of radius where 
A = 0, our analysis can be extended to other locations. 



authors take the radial axis to separate collapsing and 
expanding kinematics, here we take Eu m locally as a de- 
formed radial axis. 



B. Hypotheses and dynamical analysis 

Since part of our analysis is based on the E — Eu m 
diagram, it is useful to clarify the constraints introduced 
by the set of hypotheses we propose. 



1. Regular density distribution 

A regular density distribution is motivated by standard 
cosmological models [TTI . l29l . e.g.]. In the weak energy 
condition, the density remains positive so the mass profile 
is initially always monotonously increasing, thus Eu m , 
from Eq. 1121 is initially always monotonously decreasing, 



dM 



> 



dEu 



< 0. 



(27) 



OR ~ OR 
The regularity implies finiteness of the mass and nonzero 
values for the density at the centre. This constraints 
their logarithmic slope in the following manner: suppose 
a value — e for the slope of the density in the centre (p cx 
r _e ), then the mass shall behave accordingly as r 3-e . 
Finiteness of the mass implies then e < 3 and no vacuum 
in the centre implies e > 0, from the density. 

2. Initial Hubble-type flow 

This simplifies the initial velocity profile into one that 
only admits outgoing radial velocities (positive f) , in the 
fashion of expanding initial conditions in a Hubble flow, 
although less restrictive. As a consequence of this and 
the previous condition, the profiles in the centre always 
respect, in initial conditions, the hierarchy E < En m , 
which is crucial for the emergence of a bound core. In 
this case 



E 



(3M)s A3 

, 2AI 
~ ~R 



R 



as e < 3, 



R^O 

A 2 
3* ' 



R z 



0. 



R^0 



since r ~ R — >• so the f 2 and jR 2 both tend to 
R^> 

zero as R 2 and are thus dominated by the — ^ term for 
e > 0. Thus around the centre, 

E 2 (M 

Elim ^> _^ Q 3§i? V A 
9 

r ir 3 > i 



jz-> o 3aA3 

=> E < Ei im , for e > 0, E Um < 0. 
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In the peculiar case of a constant central density (e = 0), 
we have M ~ ^Lp i^ 3 . f ^ dnroR — H C R so 

R -+ A ->■ 

£ = {Hi - f p - A ) R-2 = - f (p + p A )) i? 2 - In 
that case, the Hubble- type flow needs to remain moderate 
in the centre to respect the constraint 

d R r < 4tt (pi (2 PA y + | (p + p A )J . 

In the rest of the paper, we assume the conditions for 
E < En rn in the centre are met. 

3. Asymptotic spatial cosmological behaviour 

If we restrict our explorations to asymptotically cos- 
mological (FLRW) solutions, this implies that at radial 
infinity the mass and velocity initial profiles, constraining 
the energy and Eu m profiles for all time, shall obey 

4tt o . , 3M 
M — > —p b R 3 with - — — — > Pb = Pb (t) 

R^oo J 4?m A^oo 

^A;™ — > - (4ir Pb )% A^R 2 , 
R -> oo 

r 4 (A) — >■ HiR^E — >■ -AA 2 . (28) 
A — > oo A — > oo 

We note that the value of the curvature A" of the 
asymptotic FLRW solution compared with the equiva- 

2 j 

lent (47rpf,) 3 As FLRW critical curvature will determine, 
together with the central constraint E < Eu m , the oc- 
curence of, at least, one intersection of E and Eu m of the 
E' > EL kind, not inducive of shell crossing (see Sec. 

MM- 

Definition 2. Supposing there exists n G N shells veri- 
fying equation hl5}) , we order them by initial radius and 
denote them A«, i € 

• R*out = R*n the outermost intersections E = Eu m of 
the initial profiles 

• R+in = R*i the innermost initial intersections E = 
Eum of the initial profiles. 

4- Local mass conservation and Lagrangian frame 

Since in our system, the cosmological constant is inert 
by definition and dust purely interacts gravitationally, 
we assume, as in [l5| . that the rest mass of each crossing 
infinitesimal shell is conserved. The shell crossing event 
can thus be viewed as an infinitesimal exchange of the 
relative positions and integrated masses while each shell 
conserves its own velocity. 

As shell masses M and energies E are conserved be- 
tween shell crossing events, Eq. (jSJ will govern the mo- 
tion of individual shells. Keeping initial R = r(0,R) 



Urn 




R 



Figure 4. Effect of an ingoing, infinitesimal test shell-crossing 
on the energy and critical energy profiles, around the local 
initial configuration for the overcoming of Eu m by E. The 
initial intersection shell becomes bound on such perturbations 
and the local intersection shell shifts outwards in radius. 



as Lagrangian labels, we can follow the dynamics of the 
shells using a the simple prescription obtained above 
without the need to reorder the radial labels as would 
require a metric extension. Instead, we keep the initial 
labels all throughout and follow each shell's evolution 
using Eq. (J5J and the shell crossing prescription of Sec. 
IIIIB41 as e.g. in Ref [HH^. 



C. Local effects of shell crossing on trapped matter 

shells 

In this section, we will detail how a test crossing shell 
affects locally the values of E and Eu m around trapped 
matter shells. 

Since each shell conserves its infinitesimal mass, the 
local effect of an elementary crossing of a system's shell 
by a test, neighbouring, shell will just exchange their 
non-local mass in the exchange of their positionso- As a 
consequence, their values of E and En m will also change. 
The change of E, in Eq. ([?]), is allowed by the shell cross- 
ing event. 

A shell crossed at some r x by an infinitesimal mass SM 
(5M > for inward crossing, < for outward crossing) 
will see its values shifted as follow (the reciprocal is true 



In this process the other shells of the system, not involved, will 
remain unaffected and conserve their masses. 



R 



Figure 5. Effect of an outgoing, infinitesimal shell-crossing on 
the energy and critical energy profiles, around the local initial 
configuration for the undercoming of Eu m by E. The initial 
intersection shell becomes unbound on such perturbations and 
the local intersection shell shifts outwards in radius. 



for the crossing shell with— SM) 



Elim+8 —Ei 



lira 



2SM 
3~TT 



-Ei 



Urn- • 



(29) 
(30) 



Thus, for an inward (resp. outward) crossing, both E 
and Eu m will decrease (resp. increase). Their relative 
separation, crucial around intersections, will follow 



SA ~26M 



ru 



(31) 



which generalises the conditions from @, 0| (see Result 
[5]). The sign of this shift is determined by the initial 
position r(to,R) = R = r.i of shells with respect to their 

'Km • 

Bound shells Cctn never cross their respective Tn 7n and 
shells with E = Eu m reach their ru m at infinity in time. 
Thus crossing events involving one bound shell, satisfy 

(-p ^-^ < 0. However, once escaping shells go be- 
yond their respective ru m , they experience the opposite 
relative effect on their SA. Thus, it is possible to have 
a crossing of two escaping shells beyond their respective 
rn m that produce shifts in the opposite direction. How- 
ever, once beyond their rn m , even drastic changes cannot 
put shells on closed orbits linked with the centre as they 
would correspond to points on the outer side of the effec- 
tive potential (Fig. QJ,). Since intersections E = Eu m 
take place in the neighbourhood of bound shells (those 
with E under their En m ) we can restrict ourselves to 
consider local shell crossing in r x < ru m . 

To first order, for inward-going crossing shells, we have 
SA < 0, as illustrated on Figs. [4] and [2 while outward- 
going shells have SA > 0, see Figs. [S] and [3] As a 



E 



5M crossing 
For r x < r Um 



r lim 



R 



Figure 6. Effect of an outgoing, infinitesimal shell-crossing on 
the energy and critical energy profiles, around the local initial 
configuration for the overcoming of Eu m by E. The initial 
intersection shell becomes unbound on such perturbations and 
the local intersection shell shifts inwards in radius. 



consequence, the limit shell defined by the intersection 
shifts forward (resp. backward) for the two cases of local 
configurations. The resulting cases are overcoming in- 
ward crossings and undercoming outward crossings (resp. 
overcoming outward crossings and undercoming inward 
crossings) and are illustrated on Figs. H] and [5] (resp. [B] 
and El). 

To simplify the qualitative study of the system, we will 
first consider a prescription where both M and E are con- 
served, in Sees. MI D ll and llll D 21 We will then drop this 
assumption and include the evolution of trapped matter 
shells' neighbourhoods, building from infinitesimal shell 
crossing as described below in Sec. IIII D 31 to ascertain 
the qualitative evolution of the system, in Sec. IIII D 41 



D. Global effect of shell crossing on limit trapped 
matter shells 



1. Simplest model with shell crossing 

In order to study the simplest set of initial conditions 
where shell crossing occurs, given the constraints of sec. 
IIII Bl from Result [51 we shall consider a model with a 
single undercoming configuration. The topological con- 
straints coming from the two dimensional E vs. R dia- 
gram together with the choice of an open background 



For example in the centre [E < Ei, 



See Fig. 
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Figure 7. Effect of an ingoing, infinitesimal test shell-crossing 
on the energy and critical energy profiles, around the local 
initial configuration for the undercoming of Eu m by E. The 
initial intersection shell becomes bound on such perturbations 
and the intersection shell shifts inwards in radius. 

at infinit\F^l leads to initial conditions for E and Eu m 
with three intersections (see Fig. [5]), the middle one ver- 
ifying Result [2] We thus have a model with R+i = i?*j n , 
R±2 and i?*3 = R+ ut defined in its initial conditions. 
We can now consider the inner system, also called the 
system, to be circumscribed by R Mu t- Unbound shells 
inside this system are in position to escape it and hence 
define a remarkable shell outside the system: 

Remark 1. The inner or non-bound shells of initial 
conditions in E and En m induce a few remarkable fea- 
tures defined as follows: 

• We will consider all shells inside R± ut as the initial 
inner system. 

• We will denote by E max the maximum value of 
non-bound E in the set of shells inside R* ut or out- 
side of it but with horizontal tangent, i.e. E max = 
max{£ : ((£' = 0) V (0 < R < R* out )) A (E > E Um )} 

• Rmax is the largest value for which E = E max , i.e. 
Rmax = max{i? : E (R) = E max } 

• Rf r ee is the furthest shell outside i?* ou t with in- 
creasing E = E max , when it exists, i.e. Rfree = 
max{i? : (R > R* out ) A (E = E max ) A (E'(R) > 0)} 

• We will note Ef ree the value of E, when it exists, as 

Efree = E (t = to, Rfree)- 

With the above definitions, we will now examine the 
effects of shell crossing on trapped matter shells. 




Figure 8. Open background with arbitrary central mass distri- 
bution and a single local undercoming intersection. It always 
gives protected inner shells as well as unmodified cosmological 
expansion, when keeping integrability despite shell-crossing. 
Shell crossing entails no fundamental modification. 

2. Limit trapped matter shells in the integrable dynamical 
system 

In a model where both M and E are conserved through 
shell crossings, we can extend the analysis of H, as each 
shell's dynamics remains integrable and is governed by 
the Lagrangian Eq. In this case, the qualitative 

dynamical behaviour of the system is entirely determined 
from the shape of its initial conditions in a E and Eu m 
vs. R diagram. 

As we will see in Sec. IIIID41 when including the full 
effects of shell crossing on E and En m , the key properties 
of trapped matter shells will be obtained in the limit 
t — > oo. Since in this section, M and E are assumed to 
be conserved with time, all the properties deduced here 
will remain unchanged in that limit. We will therefore 
express our results in the limit t — > oo, using definitions 
which evolve from Rem. [T]and are detailed in Appendix 

m 

From Fig. [51 we can see that all the bound shells 
will remain under r UmMut = r lim (R±out), while all un- 
bound shell of the inner system will escape ilF*l. Thus, 
considering that bound shells will eventually turnaround 
and orbit back and forth between the centre and their 
turnaround radius, we find that all shells inside r (t, R) = 
rum-font will be crossed from both sides (interior and exte- 
rior j^j. Only the shell i?* otl t will remain uncrossed from 
outside shells. This leads to the following definition: 



This means E — > — kp lrw R with fcFiflW < and 
R — > oo 
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The final fate of each shell will always remain on horizontal lines 
and their gravitational nature, whether bound or unbound will 
also remain the same throughout their history. 
This is indicated on Fig. |H]by horizontal arrows. 



Eum — ► ~ (47rp(j) 3 A3 ,R 2 15 This includes the R+i n and R+2 shells, locally considered trapped 

R — > oo matter shells. 
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Definition 3. The outer limit trapped matter shell 
Rt*outoo verifies Def. [7] in the limit t — > oo, in addi- 
tion to being the outermost such shell which locally is not 
shell-crossing inducive^l i.e. 



Ri 



Note, from Def. [21 that R± ut' verifies Def. [31 defines 
Rt+out, if E' > E' Um ; and verifies R tMu t = Rt*out 

oo in the 

limit t — > oo. 

Remark 2. In K-LTB with asymptotic cosmological 
evolution (FLRW at radial infinity) and initial Hubble- 
like flow (outwards going) for which shell crossing 
occurs, the outer limit trapped matter shell is a surface 
S with the following properties: 

• The matter exterior to S follows trapped geodesies, 
remaining in that exterior. 

• The matter inside S can expand and collapse protected 
from the crossing of outside shells. 

• S is the shell with largest R for which the energy 
E intersects the critical energy En m , from bound to 
unbound shells. 



The condition for existence of Rt*outoa, follows from 
the properties of Ruout, so we obtain the following result: 

Result 3. Sufficient conditions for the existence of an 
outer limit trapped matter shell are: 

• The FLRW curvature of the background k FLRW < 

2 x 

(47rpb) 3 As , or 

• Rt*out exists, or 

• The local configuration around R+out is such that 

Proof. k FLRW < (4irp b ) 3 As => E (R ->• oo) > 
Eum (R — > oo) so the last intersection E (R) = Eu m (R) 



is such that E' > 



pi 



from a corollary to Bolzano- 



Weierstrass theorem and the Definition [3] is verified. □ 

We show, in Fig. [51 a diagram with data such that the 
inner limit trapped matter shell Rt*out is at R M ut- The 
exterior of the system will include all the unbound shells 
escaping to infinity. However, the dynamics from Eq. (J5J) 
allows us to study under what conditions unbound sys- 
tem's shells will never cross shells located in the exterior 
of the svstenFl. Take two different shells R\ < R 2 , even- 



and 



tually crossing each other at a given radiui 
with the outer shell more open than the inner shell (i.e. 



E 1 < E 2 for Mi < M 2 ): 

2 A 2 2Mi 
Ei =Vi - -r x , with 



Ei <E 2l 



{R : maxji?^} A (£" (t -> oo) > E' Hm (t -> oo))} 



(32) 



16 Recall that R+oois defined by solutions in initial R of Eq. 1151 
taken at t — ¥ oo. 

17 Their escape velocity at infinity should never exceed that of ex- 
terior shells. 

18 This radius is allowed to tend to radial infinity. 



A 2 2M 2 

r 



E 2 =v 2 — — f! and 



Mi <M 2 , 



. 9 * 2AM 
=> Av 2 =AE H > and 



Av 2 - AE > 

r x — > oo 



Vt, v\ >v\. 



(33) 
(34) 

J 

(35) 



Thence shells with E 1 < E 2 and Mi < M 2 will al- 
ways remain in the same radial order and the shell with 
E = E max will then escape all other system's shells. 
It then appears that, when Rj ree exists, all shells with 
E > Ef ree will never be crossed by any shell inside i?/ ree . 
The counterpart to Def. [3] can thus be formulated by 
defining first E rnax . In turn, the condition for the exis- 



tence of E„ 



is that E > E, 



in the limit t 



Thus, Rem. [T]can be adapted here ai 



m 



Definition 4. Suppose E n 
the initial conditions, as 



exists and is defined, in 



E n 



max 



{^l(^oo): ((^ = 0) 
V (0 < R < R+outoo)) A (E > Eum) (*_►«>)} ■ (36) 
Then, inner limit trapped matter shells are defined as the 

locUS Rfreeoo Such that 

Rfreeoa = max {i? : (R > R+outoo) A (E = Emaxoo)^^ 

A(E'{t^oo,R)>0)}. (37) 

Remark 3. In K-LTB with asymptotic cosmological 
evolution (FLRW at radial infinity) and initial Rubble- 
like flow (outwards going) for which shell crossing occurs 
(and E maxoo is defined), the inner limit trapped matter 
shell is a surface S with the following properties: 

• The matter interior to S follows trapped geodesies which 
remain in that interior. 

• The matter exterior to S expands, protected from the 
crossing of inside shells. 

• S is the shell outside the system (defined with R+outoo) 
with energy equal to that of the highest E of non-bound 
shells, and starting inside of the system, or outside it but 
with horizontal tangent. 

The conditions for existence of Rfreeoo combine the 
existence of E maxoo with constraints on the background: 



19 The following can be formulated also in terms of gauge invariant 
Lie derivatives, expansion and shear, as seen in appendix [B] 
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Result 4. Sufficient conditions for the existence of an 
inner limit trapped matter shell are (a) the existence of 
Emaxoo o,nd (b) the existence of Ef reeoo : 

(a) • There exist initially a non-bound, system shell, 
or a non-bound shell with horizontal tangent: 3R : 
(0 < R < R* out V E' = 0) A E (R) > E Hm (R), or 

• R*outoa exists, or 

• There exist at least one i?« 

(b) • E maxoc < E (R -> oo), or 

• Eli? : R > R+outoo 
E'{R) > 



A E (t ->• oo, R) = E maxoo A 



Proof, (a) 



1 



If 



have R such that 



(0 < R < R* out V E' = 0) A E (R) > E Hm (R), then, 
either it is a maximum so E max exists and, by time 
evolution of its neighbourhood, E rnaxoo exists, or, by 
continuity, in the case when it is not a shell with 
E' = (local maximum), there is a shell with larger 
E which satisfies Rem. Q] for E max and thus one in its 
neighbourhood satisfying Def. [5] for E maxoo . 
2/ If R+outoc exists, it is not bound at time infinity and 
is inside the system, therefore, even if it is the only 
unbound system shell, it can at least define E maxoo . 
3/ If there is only one i?*, then it is R± ou t by Def. [2] 
We are then in the case 2/ above as this guarantees the 
existence of R+outoo- 

(b) Since: (i) E maxoo is, by definition, the largest value 
of E reached at time infinity by inner or outer local max- 
ima shells, 

(ii) uncrossed outer shells have their E conserved, 

(iii) asymptotic cosmological conditions render E 
monotonous near infinity, 

(iv) evolution of the inner shell R m axoo follows Eqs. 1351 

(v) the energy profile is continuous, 

therefore, E+ outoo < E maxoo and by continu- 
ity, since E maxoo < E(R—>oo), exterior shells 
will obey E e [-E*o«too, E (t -> oo, R -> oo)[ 2 
[Emaxoo, E (t — > oo, R — > oo)[, thus there exist at least 
one shell at time infinity with E = E maxoo . 
Moreover, for the outermost exterior shell R xm axoo = 
max{i? : R > R W utoo, E (R) = E maxoo } with E = 
E m axoo, since Emaxoo < E (R -> oo), by continuity, all 
shells outside of it will verify E > £ maxco . Therefore 

E' (Rxmaxoo) > and Rxmaxoo = Rfreeoo *S fulfilling Def. 

Ill □ 

We show a diagram in Fig. \E\ where we indicate 
the outer limit trapped matter shell for which Rf re e = 
Rfreeoo, in the model where both E and M are conserved 
between shell crossings. We thus have found, for that 
model, that extending the analysis of [Bj in the context 
of shell crossing leads to the emergence of two remark- 
able shells: an inner limit trapped matter shell and an 
outer limit trapped matter shell. From their definitions [3] 
andHJ we can deduce other properties depending on the 
background cosmological model, namely: 

From Result [3J any background with E > Eu m will 
admit an outer limit trapped matter shell. This includes 
some closed models and all flat and open models. 



From Result |H and under the assumptions of this sec- 
tion, any closed background in our models cannot foster 
an inner limit trapped matter shell as the finite value of 
Emaxoo is always larger than its energy at radial infin- 
ity. Conversely, open models always have an inner limit 
trapped matter shell (see the example of section HV| and 
only flat models with moderate enough energy fluctua- 
tions (i.e. for which E max < = E (R —J- oo)) can allow 
the existence of an inner limit trapped matter shell. In 
summary: 

Summary 1. Consider a A-LTB spacetime with asymp- 
totic cosmological evolution (FLRW at radial infinity) 
and initial Rubble-like flow (outwards going) for which 
shell crossing occurs. Then: 

1. The global limit trapped matter shell found in the 
no-shell crossing A-LTB examples of is split, if 
shell crossing occurs, into at most, two global shells, 
namely an inner limit trapped matter shell and an 
outer limit trapped matter shell. 

2. For open or flat expanding spacetimes, 

(a) there exists always an outer limit trapped mat- 
ter shell at Ruoutoo- 

(b) The inner limit trapped matter shell exists in 
flat backgrounds for sufficiently small initial 
velocities inside the system limited by R+outoo ■ 

3. For closed spacetimes, outer limit trapped matter 

2 1 

shells are present if k FLRW < (Anpb) 3 A 3 and inner 
limit trapped matter shells cannot be defined if shell 
crossing occurs, with our definitions. 

4. In the K-CDM examples of global limit trapped mat- 
ter shell found in J^j, inner and outer limit trapped 
matter shells reduce to one single surface. 

Proof. 1 : Direct from Definitions [3] and E] and Result [5] 
which leads to shell crossings at some i?* . 
2(a): From Result [3] 

2(b): Direct from Result [U as open and flat expanding 
spacetimes admit E (R — > oo) > 0. Some flat spacetimes 
can exhibit Ef reeoo > while their E — > 0. For 

R — > oo 

those cases, Definition @] is never verified. 
3: Using Result [31 for closed spacetimes, E — > — 

R -> oo 

oo <C Ef reeoo , so from Result 0] , Definition U] is never 
verified. 

4: Applying Definitions [3J and @] to configurations where 
there is only one intersection i?* = R*i = R+out verifying 
E' > E', im , no shell crossing occurs. Thus all E values 
are constant over time so R M ut = Rt*outoc , and given the 
open background, E free = S* out , so Rfreeoo = R*out = 

Rt*outoo • 

In this section, we have assumed that E and M were con- 
served through shell crossings. In the next section, we 
drop this assumption and investigate whether our previ- 
ous results remain true. 
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3. Global effect of shell crossing 

Since the sign of A = E — En m determines the binding 
property of the system, it is useful to give the final values 
of E and Eu m for each shell, labeled i, in terms of the 
initial Ri and Mi, reaching areal radius r after crossing 
shells, with 



E -> -kpLRw.R 1 



M(r(R,t),t) = Mi + / dM, 



dM n 



M, + AM, 



where the index in refers to inward crossing, out to out- 
ward crossing, j to the shells crossing shell i. 

Using definition (TT2"|) and integrating Eq. (|2TJ)) over all 
crossing shells, we get 



Eu m (r) =Ei im (Ri 



E{r) =E{Ri) - 8tt 



AM, 

M, 



3AL 



dri 



dr 



j.out 



Him 

(38) 
(39) 



where r X i is a crossing radius. Because of their qualita- 
tively simple shell crossing histories, we can look at the 
changes for three peculiar shells, singled out on Fig. [5] 
the innermost limit shell, the outermost limit shell and 
the maximum E shell initially lying in the interior of the 
outermost limit shell. 

The innermost limit shell will only be crossed by more 
bound shells exterior to it, so AMi > and 



E(r(R*i)) =E(R irl ) - 8ir J dr hm 



P( r j) r 3 



Since 



1 / AMi 



2 AMi 



3 V Mi ) V 3 M i 
AMi 



1 



Ah 



> 



2 AM 1 

< 3 Mi 



and 



< 



< 



r x i(fj) max[r x i(rj)] 



(40) 

(41) 
(42) 

(43) 



as the innermost limit shell becomes a bound shell, we 
get that 



A[£-£ H m]i <2AMi 



1 



rn m (R*i) max[r x i(r\,-)] 



< 0. 

(44) 



Thus, the innermost limit shell will globally shift out- 
wards, following the qualitative analysis of Fig. [3J 




Figure 9. Illustration, on an open background with arbitrary 
central mass distribution, of the effect of shell crossing on the 
inner global limit shell previously defined as the outermost 
local limit shell. The time variation of the locus of the outer- 
most local limit shell leads to defining it as the time infinity 
outermost limit shell: this is shown on the extrapolated time 
infinity energy profiles and linked to the initial energy profile 
by a connecting curve. The global inner limit shell is then 
just shifted inwards, compared with the integrable analysis. 



In turn, the outermost limit shell will be only crossed 
by all the unbound shells interior to it, so 



E(r(R ic0ut )) =E(R ilOUt ) + 8ir J dr^ out 



fx out i r j) 



*out )+2- 



AM out 



' xout 

)(R 

*out ) 



(45) 



where AM out is the positive mass loss of the outermost 
limit shell and (r x ) is a reduced crossing radius. Note 
that, by construction, M out > M(R i , out ,t — > oo) — 
M out ~ AM out > 0. 

Now, supposing the density distribution remains finite 
we can decompose the crossing of the outermost limit 
shell by all escaping inner shells into a series of infinites- 
imal shell crossings. Thus following Eq. (|3"Tj) we get 



d[E- E Um ] = -8-Kdrj.outp (r^ out ) r 2 jout 
1 1 



x out y j.out ) 

As all shells cross outwards^! and 

1 1 1 



> 



> 



V xout 

)(R 

*out ) Tlirn \ \> X out) ) ^lirn *out ) 



(46) 



(47) 



20 Note that R+out starts as a marginally bound shell well inside 
its limit radius. 
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freeze 
freeoo 



R > R freeoo 
never crossed 
by inner shells 



outermost R* ^Jl+out 
Elim -> - (iKPbX^R* 

Figure 10. Illustration, on an open background with arbi- 
trary central mass distribution, of the effect of shell crossing 
on the outer global limit shell previously defined as the outer 
shell with same energy function as the inner shells' maximum. 
The time variations of the inner shells' maximum energy func- 
tion from shell crossings lead to defining it as the outer shell 
with same energy as the time infinity inner shells' maximum 
energy function E: this is shown with the highest of extrapo- 
lated time infinity Es of inner shells peaks. The global outer 
limit shell is then just shifted inwards, compared with the 
integrable analysis. 




Figure 11. NFW with background En m and an example of 
E profile given by Eq. 1)501540 . setting R = 3000, arbitrarily 



setting x\ — 



X2 = 



X3 = 



ei = 10" 



= 2 and e = e 



so that the figure is proportioned. 



then, in this case, we have 
A[E- E Hm ] out = 
2AM out 



)(R 

*out ) 'Hm (R*out) 



> 0. (48) 



Thus, the outermost limit shell will shift relatively in- 
wards, following the qualitative analysis of Fig. [6l 

Finally, the maximum E shell initially inside R+ ou t , or 
with horizontal tangent (its initial radius is Rmax)-, will 
be only crossed inwards by all the shells starting with 
radii above it and having an E below E(R max ,t — > oo), 
at the moment of crossing. This shell will then follow 



A[E- E Um ] max < 2AM„ 



max[r XmM (rj)] 



< 0, (49) 



similarly as for the innermost limit shell. 

We summarize the main result of this section as: 

Result 5. Consider a K-LTB spacetime where shell 
crossing exists. Then the metric and extrinsic curvature 
are discontinuous and the discontinuity in E is given by 
Furthermore, at R M ut, A [E — E lim \ . > and, 



O-t Rmax, A [E — Eii m ] 7 



< 0. 



4- Qualitative analysis of limit trapped matter shells 

In this section, we argue that the results contained in 
Summary [T] remain true for the case where M and E are 
not conserved through shell crossing. 



We discussed the behaviour of the outermost limit shell 
R-kout and of the outward escaping highest energy in- 
ner shell R m ax in Sec. MI D 31 As those determine the 
two separating shells Rt+outoo and R freeoo studied above, 
their modifications by shell crossing will indicate that the 
effective limit shells are just displaced but obey the same 
general properties. We illustrate this on the open back- 
ground example (Fig. [8}, for which we separated the 
study of each limit shell. 

In FigEl we represent the construction of the outer 
trapper matter shell, using the qualitative evolution of 
R+out and its neighbouring shells from Eq. S5] 

In Fig llOl using the qualitative evolution of E max and 
its neighbouring shells from Eq. 1491 we represent the 
construction of the inner trapper matter shell for open 
initial conditions. The subsequent modifications proceed 
from those qualitative changes and do not modify the 
formulations of the results from their counterparts in the 
model where both E and M are conserved between shell 
crossings. 

In the case where E and M are not conserved, the 
effect of shell crossing on i?* oll t given by Eq. |48|) implies 
only that Rt±outoo < Ruout without qualitative changes 
and the Definition [3] is verified. In turn, the effect of 
shell crossing on Rf ree depends on the effect on Ef ree 
from Eq. (|49p and by the monotonous increase of E near 
infinity and only implies that Rf re eoo < Rfree- 

Therefore, the findings of Sec. MI D 21 extending the 
analysis of [5| in the context of shell crossing, are only 
quantitatively modified as full shell crossing effects only 
displace inwards both the inner and outer limit trapped 
matter shells: the initial outermost intersection of E and 
Eum gets unbound and the system at infinity gets con- 
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Figure 12. NFW with background Ey lrn and an example of E 
profile given by Eq. (|50H54[) in log (— -E)-log (-R) scale. 



sequently reduced in Lagrangian initial radius. In turn, 
the maximum energy of the inner regions gets lowered, 
so the inner limit trapped matter surface is also drawn 
inwards. This displacement modifies only marginally the 
conclusions obtained in Sec. HHP 21 namely: (i) the split- 
ting of the local trapped matter shell is maintained when 
those shells exist, (ii) Open, flat and closed models with 

existing Rt*out (with k FLRVV < (47rpf,) 3 As) all retain an 



Rt+outoo and (iii) the modification of the maximum en- 
ergy of the inner regions allows just more asymptotic 
cosmological flat models to keep their inner limit trapped 
matter shell, if the shift from E max tends to E maxoo < 0. 

Therefore, from the sufficient conditions for inner and 
outer limit trapped matter shells (Res. [f] and [3]), the 
results contained in Summary [T] remain true in the case 
where M and E are not conserved through shell crossing. 



IV. EXAMPLES: NFW PROFILES WITH ONE 
UNDERCOMING INTERSECTION 

In we studied examples of trapped matter shells 
using cosmological models with a Navarro, Frenk and 
White (NFW) density profile [3l[ and a simple parabolic 
E profile. Here we adapt those profiles in order to present 
one intersection with Eu m of the undercoming type as in 
the local configuration of Fig. [3J and thus ensure, at the 
local level, the appearance of shell crossing. 

To do so, we use a fourth order polynomial in the 
canonical Lagrange form, to provide for the behaviour 
in the intersecting region, that we cut off with an expo- 
nential so that an open FLRW term dominates at infinity. 
We chose the profile to be at the origin and at a char- 
acteristic radius Rq set near the last possible intersection 
point with the NFW E Um given by Eq. (4.21) of 0|, 

i?o < R-i,Eu m (R-i) = —1, 

so as to secure the crossings in the physical region. The 
remaining three points of interpolation are chosen to be 
set alternately below and above the En m curve inside the 
region set by and R . The form is set by 



E(R) 



mi 



x x — X2 x — X 3 x — 1 
X\ X\ — X 2 X\ — X 3 Xi — 1 



m 2 



x x — X\ x — X 3 x — 1 
x 2 x 2 —X1X2 — x 3 x 2 - 1 



m 3 



X X — X\ x — x 2 x — 1 
£3 x 3 ~ xi x 3 - x 2 x 3 - 1 



eix 



knoR 



(X 2 - 1)60 + 1' 



(50) 



where x = R/Ro, we have denoted the three intermediate 
points as x\, x 2 and x 3 , the values of the polynomial at 
those points by mi, m 2 and m 3 , t\ is a small constant 
making sure we have the freedom to fit E{Rq) = E(0) = 
where the polynomial itself is built to vanish, eo is a 



small constant making sure the polynomial dominates in 
the interesting range but allowing the curvature at radial 
infinity to be set by a Friedmann-type . The form (|50[) 
automatically vanishes at 0. We chose the polynomial 
values such that at those points, E is alternately below, 
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above and again below Eu m , the last one making sure it remains above -1: 

E(R ) =0 = eie- 1 - k^Rfo, => e = T -% 2 -, (51) 



4 

E{R X ) =gE lim {R 1 ) = ( mi + e lXl ) - k^R 2 ^ 



{x\-l)e + l 

5 



4 



(mi + eixi) e 3:1 - ^ \- — — => mi =gE Um {R 1 )e Xl 



r.3^^1 



^" U-d-^A -'J- (52) 



E\±i2) — — (m.2 + eiX2)e — ^oo^o- 



(x 2 2 - 1) e + 1 



(m 2 + eiz 2 e 3:2 + * m 2 = — -e X2 



a^e* 2 



S(i? 3 ) =E Um (R 3 ) - (E lim (R 3 ) + 1) (1 - e) 

= (E Um (R 3 ) + 1) e - 1 = (m 3 + eiafe) e~ X3 - fc^i?, 



„4 



1)60 + 1 



(m 3 + eia;3)e K3 + - - §r - ^>m3=[(75 Hm (i?3) + l)e-l]e' T 



(54) 



We illustrate this in Figs. (fTTj) and (fT2)) after choosing 
the values of the density profile identical to those in Q. 
The cut off leaves the region under Rq almost unaffected 
by the Friedmann term, so the value of is not very 
relevant there but we set it to -1. 

We summarize the result of this section as: 

Result 6. For NFW density and initial data given by 
the A-LTB spacetime has three shells i?* such 
that E\n ti = i?Hm|i? 4 - Furthermore, for this data, there 
is shell crossing and R±out — R*outoo < and Rf re e — 

reeoo ^ 0- 



V. CONCLUSIONS 

We have studied the effects of shell crossing on the 
existence of trapped matter shells in A-LTB spacetimes. 
In particular, we have considered initial conditions such 
that: (i) our models approach a FLRW solution at radial 
infinity and have an initial outgoing Hubble-type flow 
(ii) the shell crossing of dust remains pressureless and 
the mass of infinitely thin shells remains finite. 



We have shown that the local trapped matter shells 
discussed in Ref. [f| split in two shells: one outer limit 
trapped matter shell and one inner limit trapped matter 
shell. 

We have established sufficient conditions for the exis- 
tence of such shells in A-LTB spacetimes, in terms of ini- 
tial data for which shell crossing occurs. Furthermore, we 
have derived a number of properties for those shells using 
a qualitative approach inspired in newtonian-like frame- 
works of cosmological kinematical models, as in [24l. |25|. 

We have also studied the role of shear in these set- 
tings and concluded, as in [l^|, that shear favours the 
emergence of trapped matter shells. 

Finally, we have given concrete examples where shell 
crossing occurs and the inner and outer limit trapped 
matter shells emerge, using NFW data. 

As potential applications of our models we note that (i) 
Due to mass conservation and integrability in the absence 
of shell crossing at the boundary, the background asymp- 
totic conditions remain FLRW over all time. Therefore, 
this gives an interesting setting to study the extendabil- 
ity of Birkhoff 's theorem to cosmological expanding back- 
grounds; (ii) Extensions of this work to unsmooth distri- 



Rj 
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butions of mass should be possible and might give sup- 
port to current structure formation analyses using the 
spherical top hat collapse model, in the case of ACDM. 



Appendix B: Gauge invariant definitions for inner 
limit trapped matter shells 

We can rewrite E in terms of gauge invariant quantities 
withEqs. ©mm [Hand HT) 
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i 6 + 
~3 



1 (n 



r* \ r 

j2, 



E(R) — -—j^-Eii m — R z C n l — ■+ 



e 



r ' ~3 >C„ I - + a 



so the condition of existence for E max , that E > Eu m , 
translates into initial condition with the inequality 



< 



l)Ei 



lint 



R 2 



<0, 



or at time infinity into 



\ J / (t-K») r 

We get then R max and E max from 



1) E, 



! i III 



< 0. 



(t->oo) 



Appendix A: Time infinity definitions 



Definition 5. The inner or non-bound shells of initial 
conditions in E and Eu m , in the limit t — > 00, induce a 
few remarkable features defined as follows: 
• -R*oo = R*i (t — > 00) is the intersections number 
i between E (t — > 00) = Eu m (t — > 00) taken at time 
infinity but singled out by its radius in the initial profile 
of E; in particular we note R+outoo = R*n (t — > 00) for 
the outermost intersection and Ruoutoa = R*outoo when 
we add the condition {E 1 (t — > 00) > E' Hm (t — > 00)) 



We will note E, 



maxoo 



the maximum value taken at 
time infinity, but singled out in the initial profile, of 
non-bound E in the set of shells inside R+outao or 
outside but with initial horizontal tangent, i.e. E maxoo = 
{E, max (E (t -> 00)) A ((£' = 0) V (0 < R < iWoo)) 
A (E > E Hm )} 
• R, 



■maxoo 



is the largest value for which E = E ri 



R 

maxoo — 

max{i?, E(R) = E maxoo \ 

• Rfreeoo, if it exists, is the furthest shell outside R* utoo 
with an increasing E at E = E maxoo , i.e. Rfreeoo = 
max{i?, (R > R^outoo) A (E = E rnaxoo ) A (E'(R) > 0)} 

• We will note Ef reeoo the value of E. if it exists, such 

0,S E freeoc — E (T — 0, Rfreeoo^) • 



Rmax = max < R, = - min lc n ( ® + a 



Ei 

R3~ 



E„ 



A((E' = 0) V (0 < R < R* out ))f\ I C n + 



< 



DE-, 



I ill) 



R 2 



< 



i i i En m rn m 
-R mnn . mm <£„(— + a 1 — 



max < R- 



2 ( Elijnrii 



R-- 



Taken at t — > oo, this translates into 



((£' = 0) V (0 < R < R* out )) 

<oH. 



G 




Rmaxoo = max { R, E = max {-r (£„(— + « 



ElimTli 



A {(E' = 0) V (0 < R < R wut0 o)) 



(t— >oo) 




(t->oo) 
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Definition 6. Suppose that E maxoo defined as 



Emaxoo — max \ T 



max < r 



ElimTli 



' ¥ + a 



E, 



I 6 a. 
-£„(—+ a 



(t->oo) 



(t-S-oo) 



((£' = 0) V (0 < R < R* out oo)) 



e 



A £„ - + a < — 



< 



(t->-oo) 



((£' = 0) V (0 < i? < i?» outoo )) 



( I + a) < — 



(t-S-oo) ^ 



(Bl) 



exists. Then, inner limit trapped matter shells are de- 
fined, in the models considered with GLTB coordinates, 
as the locus Rf ree oo such that 



Rfreeoa ~ max j R, (R > R* ou tco) 

( e 



3 



En 



M 1 



A 



A (E'(t -> oo, i2) > 0) L (B2) 



Thus, Def. @]can be rewritten as 
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